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SUMMARY 

A computer program was developed to calculate numerically the speed and 
circling polars of an aircraft when the lift and drag characteristics of the 
wing airfoils are known. The planform of the wing is described by variables 
which are optimized so that the cross-country speed of the glider is maxi- 
mum for the particular type of thermal model. Two thermal models will be 
compared and it can be shown that with a greater wing area than now normally 
used the performance can be increased. 
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SYMBOLS 

half span 

section drag coefficient 

section lift coefficient 

total drag coefficient of the aircraft 

total lift coefficient of the aircraft 

induced drag coefficient 

parasite drag coefficient 

acceleration of gravity 

height 

radius of thermal 

area 

time 

net climb rate of the glider 
sink rate for the speed polar 
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sink rate for the circling polar 
vertical velocity of thermal 
airspeed 

average cross-country speed 

cross-country speed for the i-th thermal 

total weight 

object function 

variables 

weighting factors 

air density 

sweep angle 

aspect ratio 

bank angle 


INTRODUCTION 

The most important part of an airplane is a well designed wing. Only a 
few variables are necessary to describe the planform (fig. 1). The halfspan 
b and the sweep-angle <j> are fixed. The airfoils with their lift and drag 
characteristics must be prescribed. The total lift of the airplane is 
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( 1 ) 


The C is calculated either by Multhopp’s method (ref. 1) or, for swept wings, 

J-t 

by Truckenbrodt ' s method (ref. 2). The tail lift is in general so small that 
it can be neglected. From the lift calculation the induced drag (ref. 1) 

is obtained, too. The total drag is i 


C D “ C D + S. + C Do 

profile i 
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is the drag coefficient of the p-th part, with surface area S . The 
Do p P 


profile drag coefficient C 


can be calculated from the measured 


profile 


Reynolds-number-dependent c^-versus-c^ plots (fig. 2). 
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From fig. 2 (measured data), fig. 3 can be determined, and by linear interpo- 
lation, the local c. at the strip n with surface S is obtained. The air- 
* d v n 


n 

craft equivalent parasite area y Sp is assumed to be constant for a 

P P 


given aircraft. 


SPEED AND CIRCLING POLARS 

With the weight as the fifth variable (x^ = weight) , the principal per- 
formance characteristic of the glider, which is its cross-country speed, can 
be calculated. The speed polar is given by the equations 
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CROSS-COUNTRY SPEED 


The net rate of climb of the glider in the thermal (fig. 4) is 


v (r) = v , .. 

c thermal 


(r) - v (r) 
sc 


( 9 ) 


The maximum climb rate must be calculated as a function of the diameter of the 
thermal. With these two polars (speed and circling) the cross-country 
speed can easily be determined (fig. 5) , 
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(D = distance; t = time from point A to C) . The loss in height is 

H - t _ v 
1 s 


The gain in height is 


H = 


• v 
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Therefore, 
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(13) 


and with the time t^ to go from point B to C, and the total time t, re- 
spectively, 
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(15) 


the cross-country speed by using the equations (13) to (15) and (10) is 
finally 

V • v 

V = — — 

R v + v 
c s 


(16) 


THE OBJECT FUNCTION 

The program can be used with different types of thermal models. Two 
examples are given in fig. 6 and fig. 7 (ref. 3,4). As object function for 
the optimization, an average cross-country speed is defined for each thermal 
model 
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Z = Ea/V.^ (i = 1,2,3 for Carmichael) (17) 

i i (i = 1,2,.. 4 for Horstmann) 

U. = 1 and 1 j> ct. _> 0 (18) 

i 

The V ’s are the cross-country speeds from the i-th thermal and a. is its 
R i 1 

weighting function. The ou f s must be chosen by the designer and then the 

optimization is achieved for the particular distribution of thermals assumed. 


THE OPTIMIZATION PROCEDURE 

The maximum cross-country speed is calculated by a penalty optimization 
program (ref. 5,6), so that constraints like x. > a. (i = 1,2,.. 5) (a. = 

i—i i 

lower bound) or others can be observed. Using the performance polars 
from equations (1) and (2), two sub-optimization problems are 
solved: (a) determination of max v^ by using C^ and r as variables in 

equations (7) to (9), and (b) determination of max V (eq. 16) by using 

Cp as a variable through equations (5) and (6), with max v^ fixed. V R 

has to be calculated for each thermal (index i) . Then Z (eq. 17) can 1 
be calculated and optimized. 


RESULTS 

Only gliders with a span of 15m = 2b and $ = 0° were optimized. In 
diagrams 1 and 2 an optimization is shown with the profile FX-61-184 (ref. 7) 
for two thermal models with different sets of a. T s. For a convenient repre- 
sentation (this is only one possibility) , was used for the 

Carmichael thermals and a 1 = a n and a„ = a. for the Horstmann thermals. For 

11 3 4 

comparison, the performance of the D-38, an almost optimized competition 
glider in the 15m class from the Akaflieg Darmstadt (ref. 8), is shown in 
both diagrams. The curves were calculated by changing only the weight of 
the glider. The difference between the two curves gives the gain in per- 
formance over the D-38. On the average, the gain in performance is more 
than 2% for both thermal models. The Carmichael thermals deliver a higher 
cross-country speed because it is narrow and strong. This leads to a higher 
optimal wing area for Carmichael thermals compared to the Horstmann thermals 
(see table 1 and 2). In general, optimized wing area is significantly 
larger than the normally used 11m area for rigid 15m wings. 
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In diagram 3 the wing loading and the aspect ratio are plotted versus 
the variable sets of a^'s for both thermal models, using the numbers from 

tables 1 and 2. The difference in the aspect ratio, AX , is about 1 and 
therefore the wing area difference, AS, is about 1m 2 . The average wing 
area is 14.4m 2 for the Carmichael and 13.5m 2 for the Horstmann thermals. 
Clearly the influence of the shape of the wing planform is small. The 
aspect ratio is almost constant from weak to strong thermals. This is true 
for both thermal models. The wing loading is the dominant factor influencing 
the cross-country speed. The dotted line (diagram 2) which was obtained by 
changing the weight, using optimal design for all a^'s equal, shows the im- 
portance of the weight (wing loading) but a relatively small influence of 
the wing shape. The proper wing loading, again, is more important for the 
Horstmann model (greater gradient) than for the Carmichael model for 
maximum cross-country speed. 

Because of the dominating rule of the weight for the gliders, the 
curves in diagram 4 were calculated for various masses by using the 
optimal design with all OL f s equal. The optimal variables are given 

in table 3. The comparison between the 7 profiles shows that the best rigid 
profile is the FX-61-184, and the best profile with flaps is the FX-K-170 


CONCLUSION 

An increase of performance of more than 2% is possible with an increase 
in wing area. It is not necessary to design a glider for extreme (strong 
and weak) thermal conditions. The use of average weather (all a^’s equal) 

for the optimization of a glider results in an almost optimal design for 
all weather conditions. But the glider must be built as light as possible 
and should be able to carry up to 150 kg water ballast. The thermal model 
(Carmichael, Horstmann, or perhaps others) plays a minor part in the design 
but is extremely important for the proper choice of the water ballast to 
maintain maximum performance. Here it matters whether to believe in the 
Carmichael, Horstmann, or other thermal models. The program can then be 
used for a single variable optimization (the weight) to calculate the 
optimal water ballast curves for any type of weather condition. 
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Table 1: 


Optimal Variables for 15m Gliders with the Profile FX-61-184 
(Carmichael thermals, C^S = 0.04m^) 


OPTIMAL 

OPTIMAL 15m GLIDERS 

D-38 

VARIABLES 

“l =a 2 =0 

°3 =1 

a 3 =a 2 = 0. 2 
<* 3 = 0.6 

a l a 2 a 3 

a l = a 2 = ^ ’ 45 

o 3 = 0.1 

a l =a 2 = 0 * 5 
a 3 =0 


(m) 

1.31 

1.19 

1.23 

1.14 

1.19 

0.94 

x 2 (m) 

1.11 

1.07 

1.09 

1.04 

1.02 

0.753 

x 3 (m) 

0.49 

0.47 

0.51 

0.45 

0.45 

0.376 

(m) 

3.78 

4.17 

3.87 

3.62 

4.25 

4.5 

x 5 (kg) 

249 

289 

327 

360 

381 

300 

DERIVED 

VALUES 







s (m 2 ) 

14.89 

14 . 6 

14.76 

13.65 

14.2 

11 

X 

15.11 

15.4 

15.22 

16.48 

15.85 

20.45 

max C L /C D 

34.3 

34.8 

34.8 

35.7 

35.5 

37 
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Table 2: Optimal Variables for 15m Gliders^ith the Profile FX-61-184 

(Horstmann thermals, C S = 0.04m ) 


OPTIMAL 

VARIABLES 


OPTIMAL 

15m GLIDERS 


D-38 

a i =a 2 = 0 

a =a =0. 5 
3 4 

a l~ a 2~0 * 25 

a =a =0. 25 
3 4 

a r“2 =0 - 4 

a =a =0. 1 
3 4 

a l =a 2~0 * 5 

a =a, =0 
3 4 


(m) 

1.32 

1.19 

1.149 

0.999 

0.94 

x 2 (m) 

0.99 

0.884 

1.012 

0.995 

0.753 

(m) 

0.507 

0.384 

0.394 

0.389 

0.376 

x 4 (m) 

4.05 

4.48 

4.28 

3.679 

4.5 

X 5 (kg) 

254 

328 

390 

432 


DERIVED 






VALUES 






S (m 2 ) 

14.53 

13.2 

13.5 

12.64 

11 

X 

15.48 

17.04 

16.67 

17.8 

20.45 

max C L /C D 

34.6 

35.9 

36.0 

36.7 

37 
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Table 3: Optimal 15m Gliders with Different Profiles 

(Horstmann thermals, C^S = 0.04m2) 



«i a 2 a 3 a 4 

Profile 
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OPTIMAL 
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rH 

1 

vO 

l 

rH 

1 

i 

i 

VARIABLES 

vo 

1 

VO 

1 

VO 

1 

vO 

l 

¥ 

¥ 

X 

1 


►x} 

X 

X 

X 

X 

X 



Fh 


X 

F 1 -! 

Fh 

Fh 

Fh 

x (m) 

X 

0.94 

1.07 

0.98 

1.19 

1.18 

1.03 

1.077 

x 2 (m) 

0.76 

0.927 

0.919 

0.884 

0.99 

0.866 

0.911 

x 3 (m) 

0.379 

0.395 

0.37 

0.384 

0.38 

0.39 

0.37 

x^ (m) 

4.37 

4.08 

3.37 

4.48 

4.9 

4.26 

4.14 

x 5 (kg) 

290 

311 

316 

328 

348 

344 

367 

DERIVED 








VALUES 








S (m 2 ) 

10.76 

12.8 

11.52 

13.2 

14.08 

12.1 

12.74 

X 

20.9 

17.58 

19.53 

17.04 

15.98 

18.6 

17.66 

max C L /C D 

36.4 

35.7 

34.3 

35.9 

38.5 

37.6 

37.1 
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Diagram 1: Optimal performance of a 15m glider with the profile FX-61-184 

for the Carmichael thermals. The numbers on the curves are total 
weights of the gliders in kg. (C^S = 0.04m 2 = const.) 
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Diagram 2: Optimal performance of a 15m glider with the profile FX-61-184 

for the Horstmann thermals. The numbers on the curves are total 
weights for the gliders in kg.' (C* S = 0.04m 2 = const.) 
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Diagram 3: Comparison Between Carmichael and 

Horstmann Thermals. 



Diagram 4: Performance of 15m gliders using different profiles for the Horstmann 

thermals. Curves are calculated by using the optimal design ( ^ design 
point) with all a . ' s equal. The change of the total mass from point 
to point is 50 Kg^C^ S = 0.04m 2 ) 
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